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CHAPTER  I 


INTRODUCTION 

Of  the  various  processes  necessary  to  produce  lasing  from  a  chemical 
laser  device,  one  of  the  critical  processes  in  the  H2- ^2  continuous  laser 
is  the  reactant  mixing  in  the  optical  cavity.  It  is  known  that  the  chemical 
laser  performance  is  dependent  on  the  relative  rates  of  excited  state  produc¬ 
tion,  collisional  deactivation,  and  the  radiative  processes.  Consequently, 
it  is  clear  that  the  mixing  rate,  through  its  influence  on  the  production 
rate,  plays  an  important  role.  It  was  the  objective  of  the  work  discussed 
in  this  report  to  investigate  experimentally  the  mixing  process  and  its 
influence  on  laser  performance. 

Chemical  laser  systems,  currently  under  investigation,  involve  extremely 
complicated  interactions  of  chemistry,  kinetics,  optics  and  fluid  dynamics. 

Of  primary  concern  in  the  laser  chemistry  arethe  excited  states  of  species 
in  the  laser  cavities  that  provide  the  desired  radiant  energy.  For  many 
systems  the  exchange  or  transfer  of  excited  states  from  one  specie  to  another 
and  the  need  for  rapid  combustion  and  interactions  in  the  near  field  is 
essential  for  producing  the  desired  radiation.  Current  applications  involve 
the  high-speed  exhaust  of  hot  combusting  gases  that  are  expanded  through  a 
two-dimensional  contoured  nozzle  system  with  alternate  nozzles  of  dissimilar 
gases  for  vibrational  energy  transformations.  Parallel  sonic  and  supersonic 
jets  are  currently  being  used  to  inject  the  chemicals  into  the  cavity  for 
mixing  and  power  generation.  To  maximize  the  energy  from  the  laser  the 
gases  must  mix  rapidly  in  the  near  field  to  provide  rapid  molecular  and 
atomic  diffusion  of  the  gases.  Fluid  dynamic  mixing  in  the  laser  cavity  is 
the  forcing  agent  that  places  the  gaseous  species  in  a  position  for  maximum 
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and  rapid  diffusion. 

The  research  effort  described  here  was  limited  to  an  investigation 
of  the  fluid  dynamic  mixing  of  the  gases  as  they  enter  the  near  field  of  the 
laser  cavity,  to  correlate  theory,  test  data  and  photographs  of  the  mixing 
phenomena,  and  to  establish  the  effects  and  sensitivity  of  properties  on 
mixing  of  parallel  supersonic  streams.  To  this  end,  holographic  impres¬ 
sions  were  obtained  for  cold  flow  tests  and  were  used  to  reconstruct  flow 
field  views  (shadowgraphs,  schlierens  and  i nterferograms )  for  use  in  the 
evaluation  and  analysis  of  the  cavity  flow  field. 
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CHAPTER  II 

CHEMICAL  LASER  SYSTEMS 

Polanyi^  presented  many  of  the  basic  concepts  we  know  today  for  the 
vibrational -rotational  energy  transfer  of  the  chemical  laser  phenomenon. 

In  his  paper  he  observed  that  vibrational  excited  states  have  long  sponta¬ 
neous  emission  lifetimes  and  that  for  many  chemical  reactions  the  energy  of 
reaction  is  converted  efficiently  into  vibrational  states.  He  defined  the 
concepts  of  complete  and  partial  vibrational  population  inversions  and  showed 
that  both  would  occur  as  a  result  of  reactions.  The  possibility  of  an 
exchange-reaction  laser  consequently  existed  and  it  was  suggested  to  use 
low  pressure,  low  temperature  and  fast  flow  to  minimize  collisional  deacti¬ 
vation  and  thus  provide  a  continuous  chemical  laser. 

2 

Kasper  and  Pimental  later  reported  on  the  concept  of  the  vibrational 
transition  laser,  which  works  on  a  simple  chain  principle.  After  the  produc¬ 
tion  of  a  few  Chlorine  atoms  by  flash  photolysis,  the  chain  proceeds  until 
all  of  the  initially  available  H£  and  Cl 2  molecules  have  been  consumed.  The 
first  step  in  the  chain,  being  slightly  endothermic,  produces  only  ground- 
state  HC1 .  The  second  step,  however,  being  considerably  exothermic  produces 
HC1  vibrational  levels  to  about  the  sixth  level. 

The  chemical  formulation  involving  atomic  fluorine  and  ammonia,  atomic 
and  evolution  hydrogen  and  chlorine,  fluorine  and  hydrogen,  etc.,  typically 
used  in  chemical  lasers  was  not  addressed  specifically  in  the  experimental 
program  except  to  simulate  the  molecular  weight  and  specific  tests  of  the 
products  in  a  cold  flow  simulation  of  a  continuous  supersonic  diffusion 
laser  (SDL).  In  this  chemical  laser  device,  D£  is  used  in  the  combustion  of 
an  HF  laser,  and  vice  versa,  to  prevent  absorption  of  the  laser  radiation 
in  the  combustion  products.  Several  types  of  nozzle/injector  configurations 
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have  provided  good  lasing  performance,  and  were  the  basis  for  selecting 

the  various  nozzle  configurations  for  the  cold  flow  tests. 

3 

As  discussed  more  fully  by  Warren  ,  Figure  1  shows  a  simplified  sketch 
of  the  flow  downstream  of  the  nozzle/injector  system  in  a  combustion-heated 
SDL.  Wedge-shaped  nozzles  are  shown,  but  contoured  nozzles  are  also  used. 

For  clearness,  the  sketch  is  for  planar,  or  slit,  nozzle  geometries.  Because 
of  the  small  nozzle  size  (2  to  12  mm  exit  height)  and  generally  low  flow 
Reynolds  numbers,  fairly  thick  wall  boundary  layers  are  formed.  Some  systems 
have  been  operated  with  fully  viscous  nozzle  flows;  that  is,  with  a  signif- 
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leant  degradation  of  centerline  total  head. 

The  flows  exiting  from  the  individual  nozzles  and  injectors  are  confined, 
except  near  the  edges  of  the  nozzle  injector  bank,  by  the  flows  from  the 
neighboring  nozzles  and  injectors.  Thus,  the  main  flow  is  nearly  a  constant- 
area  process  composed  of  a  repetitive  series  of  similar  flows  with  short 
characteristic  dimensions  along  the  optical  axis.  Because  of  the  high  energy- 
release  rate,  the  overall  conservation  laws  require  a  significant  average 
static  pressure  rise.  This  is  achieved  through  the  generation  of  shock 
waves,  which  originate  at  the  nozzle  lips  and  initially  propagate  into  both 
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the  unreacted  fluorine  and  streams. 

Each  shock  influences  the  flows  from  several  adjacent  no;:zles  while 
reaction  and  lasing  are  occurring.  The  strength  of  the  initial  shock  waves 
and  their  effects  on  the  unreacted  stream  properties  and  flow  directions 
have  been  estimated  through  a  displacement-thickness  analysis  and  been 
found  to  provide  significant  perturbations  to  the  initial  undisturbed  proper¬ 
ties.  Such  processes  have  not  yet  been  included  in  theoretical  descriptions 

3 

of  a  lasing  flow. 

Downstream  of  the  nozzle/injector  face,  the  fluorine  and  streams 


Figure  1.  Description  of  the  flow  processes  in  a  supersonic 
diffusion  HF  chemical  laser  (Reference  1). 
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diffuse  into  each  other,  and  the  reaction  process  is  initiated  by  the  F 
and  H2  reaction.  In  most  SDL  systems,  a  considerable  excess  of  ^  (of  H2 
plus  diluent)  is  used  compared  to  stoichiometric,  and  these  streams  will 
appear  to  persist  further  downstream  than  do  the  unmixed  cores  of  the  fluorine 
streams.  The  reacting  diffusion  zones  formed  are  essentially  free  shear 
layers  generated  by  the  momentum  of  velc  'ity  differences  between  the  two 
initial  streams  and,  to  a  first  approximation,  may  be  assumed  to  be  governed 
by  the  reacting  boundary-layer  system  of  equations  (away  from  the  wave- 
interaction  regions).  The  diffusion  process  may  be  laminar,  transitional, 
or  turbulent;  and  it  is  possible  but  not  yet  established,  that  all  three 
processes  significantly  influence  the  flow  in  present  SDL  devices. 

For  a  given  geometry,  an  increase  in  the  overall  pressure  level  of  tne 
flow  can  adversely  affect  lasing  efficiency  in  several  ways.  The  diffusion 
zones  will  be  thinner,  at  least  for  the  laminar-mixing  case,  and  this  will 
delay  the  reaction  of  all  the  fluorine.  The  two-body  collisional  deactiva¬ 
tion  exchanges  will  also  be  faster.  Such  considerations  have  been  modeled 
for  a  SDL  using  constant-pressure,  flame-sheet,  and  two-vibrational-level- 
molecule  assumptions.  Scaling  laws  have  been  developed  for  both  laminar 
and  turbulent  mixing  in  the  reaction  zones,  that  quantify  the  behaviors 
involved.  Similar  results  have  also  been  developed  for  a  multi-vibrational- 
level  model.  Reynolds  scaling  can  be  applied  to  the  SD!  lasing  process, 
since  the  im;ortant  chemical  production  and  deactivation  processes  are  all 
two-body  react iors  and  since  the  mixing  processes  are  Reynolds-number  depend¬ 
ent.  The  couplea  reaction  flow  processes  will  scale  if  the  product  of  den¬ 
sity  and  length  in  the  system  is  held  constant.  Such  scaling  is  quite  gene¬ 
ral  in  that  it  governs  all  of  the  complex  processes  in  the  flow.^ 


Chapter  III 
HOLOGRAPHIC  SYSTEM 


The  holographic  and  reconstruction  systems  used  in  the  analyses  and  diag¬ 
nosis  of  the  laser  flow  fields  was  a  basic  optics  network,  including  three 
lasers,  reflectors,  beam  splitters,  collimators,  and  spatial  filters,  which 
collectively  make  up  the  systems.  Details  of  the  specific  system  is  discussed 
in  the  chapter  concerning  the  experimental  program.  Other  information  on  the 
fundamentals  of  the  shadow  method,  Schlieren  technique  and  interferometry  is 
presented  in  this  chapter  and  in  the  Appendix. 

Basically,  the  hologram  of  a  three-dimensional  scene  may  be  considered  as 
a  window  through  which  a  precise  image  of  the  scene  can  be  viewed.  The  radia¬ 
tion  field  on  the  viewer's  side  of  the  hologram  has  the  same  information  con¬ 
tent  in  either  case.  This  radiation  field,  which  is  unique  to  a  given  object, 
can  be  extremely  complex.  It  may  be  reconstructed  by  illuminating  the  hologram 

with  a  simple  wave  front  which  through  the  process  of  diffraction  is  converted 
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into  a  more  complex  one.  The  coherent  radiation  field  is  referred  to  here, 
which  originally  is  reflected  from  or  passed  through  the  object  field.  Any 
radiation  generated  by  the  object  field  must  be  excluded.  Such  radiation  rarely 
meets  coherence  requirements  and  it  simply  adds  a  bias  to  the  coherent  radia¬ 
tion. 

Holograms  can  be  made  by  many  methods  on  many  types  of  material.  The 
methods  applied  here  were  limited  to  optical  holograms  produced  on  photographic 
emulsions.  A  hologram  can  be  made  by  recording  the  intensity  of  the  sum  of 
two  or  more  waves.  One  of  the  waves  must  be  a  reproducible  wave  (the  reference 
wave).  The  other  is  the  modulated  or  object  wave  carrying  phase  and  amplitude 
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information  impressed  upon  it  by  an  interest  field  (the  object  wave).  The 
resulting  recorded  interference  pattern,  or  grating,  is  the  hologram. 

Except  for  simple  cases  a  hologram  is  a  complicated  diffraction  grating 
much  too  intricate  to  produce  by  methods  usually  employed  to  make  dif¬ 
fraction  gratings.  Their  reproduction  is  also  made  possible  through  inter¬ 
ference  of  coherent  waves,  and  subsequent  recording  of  the  interference 
pattern.  A  hologram  of  a  point  object  is  a  grating  which  when  illuminated 
by  coherent  light  will  diffract  a  significant  portion  of  that  light  to  a 
point  in  space.  The  grating  consists  of  a  series  of  concentric  bright  and 
dark  rings.  Holograms  of  more  complicated  objects  are  then  understood  as 

linear  superposition  of  many  such  diffraction  gratings  which  produce  points 

4 

in  the  image  of  the  object. 

If  a  photographic  plate  is  exposed  to  a  field  of  fringes  created  by  the 
interception  of  laser  light,  the  resulting  recordings  form  the  basis  for  holo¬ 
gram  classification.  The  developed  emulsion  contains  surfaces  of  deposited 
silver  along  the  bright  fringes.  If  the  developed  emulsion  is  returned  to  its 
original  recording  position  and  illuminated  by  only  one  of  the  original  point 
sources,  the  light  reflected  by  the  planes  in  the  emulsion  will  appear  to  an 
observer  to  emerge  from  the  original  position.  When  used  with  the  hologram, 
either  of  the  point  sources  can  be  used  to  reconstruct  an  image  from  the  other. 
A  three-dimensional  object  source  in  general  is  made  of  many  composite  points. 

The  reference  source  is  usually  a  point  source  or  collimated  source  of  radia- 

4 

tion  which  can  be  duplicated  for  the  reconstruction  process. 
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In  order  for  the  fringes  to  be  distinct,  the  object  and  reference  wave- 
front  source  must  be  fixed  relative  to  each  other  to  such  extent  that  the 
fringe  planes  will  remain  fixed  in  the  emulsion  during  recording.  The  fringe 
planes  cannot  be  allowed  to  move  more  than  about  1/10  the  separation  distance 
between  fringe  planes.  The  criterion  is  most  severe  where  motion  of  either 
source  exists  and  will  deteriorate  the  recording.  The  criterion  is  most 
easily  met  where  motion  greater  than  A  is  tolerable.  The  film  resolution  cap¬ 
ability  must  also  be  at  least  equal  the  fringe  spacing.  Again,  the  white  light 
case  is  most  severe,  requiring  A/2  resolution.  For  green  light  this  is  in  the 
order  of  5,000  lines  per  millimeter,  approaching  the  limit  of  the  best  photo¬ 
graphic  emulsions  available.  As  the  angle  between  the  object  and  reference 
beam  decreases,  so  does  the  film  resolution  requirement,  reaching  its  minimum 
for  the  in-line  case. 

Except  for  the  white  light  reflection  case  a  hologram  should  be  made  with 
its  normal  bisecting  the  angle  between  object  and  reference  radiation  to  reduce 
the  effects  of  emulsion  shrinkage.  Emulsion  shrinkage  is  largely  normal  to  the 
plane  of  the  emulsion.  If  the  above  criteria  is  met,  fringe  planes  will  lie 
normal  to  the  emulsion  and  shrinkage  will  not  change  the  fringe  spacing.  This 
effect  becomes  more  severe  for  thicker  emulsions.  When  the  emulsion  thickness 
is  increased  to  the  same  order  as  the  fringe  spacing,  the  hologram  must  be 
illuminated  by  the  reference  wave  more  precisely  at  the  same  angle  used  in 
recording.  This  allows  more  than  one  image  to  be  stored  in  a  thick  emulsion 
hologram  and  if  different  reference  beam  angles  are  used,  the  images  can  be 
reconstructed  independently.  When  a  thin  emulsion  hologram  is  illuminated  by 
reference  source  not  identically  matching  the  original  reference  source  used 
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to  make  the  hologram,  three  things  happen:  (1)  the  reconstructed  image  po¬ 
sition  changes,  (2)  the  reconstructed  image  exhibits  magnification,  and 
(3)  the  image  becomes  aberrated.  Changing  the  radius  of  curvature  of  the 
reference  wave  causes  spherical  aberration.  Changing  its  angle  of  illumina¬ 
tion  of  the  hologram  causes  distortion,  coma,  and  astigmatism.  Changing  its 
wavelength  causes  chromatic  aberration.  The  effects  of  aberrations  become  more 
serious  as  fringe  spacing  is  reduced.4 

The  utility  of  holography  ultimately  depends  upon  the  ability  of  an  in¬ 
vestigator  to  derive  data  from  the  reconstructed  image.  When  a  hologram  of  an 
extended  volume  of  many  particles  is  analyzed,  any  individual  particle  image 
is  emersed  in  light  scattered  from  out-of-focus  images  in  addition  to  optical 
noise.  Since  coherent  optical  noise  often  has  the  same  appearance  as  smaller 
particle  images,  considerable  skill  is  sometimes  required  to  distinguish  them. 

As  far  as  film  properties  are  concerned  for  holographic  purposes,  it  has 
been  shown  that  smaller  particle  sizes,  larger  particle  distances  from  the 
recording,  and  higher  resolution  of  the  reconstructed  image  require  higher 
resolution  film.  For  every  holographic  recording,  the  effective  diffraction 
pattern  diameter  is  limited  either  by  the  film  noise,  resolution,  physical 
size  or  a  combination  of  these  factors.  This  effective  diameter  and  the  dis¬ 
tance  of  the  object  from  the  hologram  define  the  recorded  solid  angle  of  light 
and,  therefore,  the  maximum  resolution  in  the  reconstructed  image.  Films  are 
available  that  minimize  the  above  factors.4 

Silver  halide  photographic  emulsion  on  glass  is  typically  used  in  holog¬ 
raphy.  The  choice  of  emulsion  is,  however,  more  critical  for  success  in  hol¬ 
ography  than  with  photography.  Choosing  a  film  with  too  low  of  resolution 
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for  a  particular  problem  in  photography  will  lead  to  a  substandard  image.  The 
same  error  in  holography  may  lead  to  a  recording  with  no  capability  to  produce 
an  image  at  all.  The  primary  factors  to  be  considered  are  resolution  and  re¬ 
quired  exposure.  Other  factors  are  concerned  with  noise  properties,  and  special 

4 

absorption  and  scattering  properties. 

The  developed  hologram  can  be  reconstructed  by  illuminating  the  hologram 
with  laser  light,  preferably  at  or  near  the  frequency  of  which  the  hologram  was 
made. 

Reconstructed  images  may  first  be  scanned  under  low  magnification  to 
provide  a  quick  examination  of  the  volume.  Portions  can  then  be  enlarged 
by  increasing  the  distance  from  the  lens  to  the  camera  and  maintaining 
focus  with  the  hologram.  Noise  from  irregularities  on  glass  surfaces  in  the 
reconstruction  system  usually  limits  useful  lens  magnification  to  about 
ten. 

The  optics  network  can  also  be  used  to  make  interferograms .  The  basic  . 
principles  of  holographic  interferometry  are  not  different  from  conventional 
interferometry.  Nevertheless,  it  is  this  field  that  has,  perhaps,  been  ex¬ 
panded  more  so  than  any  other  by  the  introduction  of  holography.  As  discussed, 
interferometry  is  the  mixing  of  two  or  more  coherent  waves  for  the  purpose  of 
making  measurements  which  are  manifested  in  the  resulting  interference  fringes. 
Holographic  interferometry  is  merely  interferometry  in  which  some  of  the  asso¬ 
ciated  waves  are  stored  holographically.  In  one  class  of  interferometry,  a 
comparison  wave  of  known  character  is  mixed  with  a  test  wave  of  unknown  charac¬ 
ter  to  determine  its  spatial  phase  variation  impressed  thereon  by  a  test  field 
or  object.  In  a  second  class  of  interferometry,  two  or  more  waves  of  unknown 
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character  are  mixed  for  the  purpose  of  determining  the  spatial  phase  variation 
of  the  two  waves.  In  any  of  these  cases  holography  provides  a  storage  medium 
for  comparison  waves.  Typically,  test  waves  which  exist  at  some  instant  in 
time  are  stored  for  late  recall.  This  provides  for  mixing  and  comparing  them 
with  comparison  and  test  waves  which  are  variable  at  will  in  shape,  wave¬ 
length,  and  mixing  geometry.  It  also  provides  more  time  for  analysis  of  the 
test  wave  which  existed,  an  important  factor  in  fluid  dynamic  cases. 

The  type  of  holography  used  in  this  experiment  was  one  of  direct  illum¬ 
ination  where  laser  light  is  passed  through  a  flow  field  test  cavity  and  is 
modulated  in  phase  by  the  gas  density  gradients.  Such  a  beam  can  be  recorded 
in  a  hologram,  and  its  analysis  can  be  carried  out  at  a  later  time.  This  pro¬ 
vides  for  a  "frozen  time"  optical  analysis  of  the  field.  Any  of  the  conven¬ 
tional  flow  visualization  methods  can  be  performed  on  the  reconstructed  wave. 
These  include  shadowgraph,  schlieren,  interferometry  and  many  other  phase  sen¬ 
sitive  methods.  From  a  single  hologram  made  at  a  given  instant  in  time,  one 
has  an  almost  endless  choice  of  viewing  schemes  which  can  be  employed  to  examine 
specific  flow  parameters. 

Both  the  pulsed  lasers  and  CW  lasers  can  be  used  in  forming  the  hologram, 
while  continuous  wave  lasers  are  used  for  the  reconstruction.  The  adaptation 
of  a  conventional  flow  visualization  system  for  holography  requires  the  use  of 
a  laser  illuminator,  and  the  addition  of  a  reference  beam.  The  hologram  com¬ 
monly  records  the  wavefront  converging  from  the  second  primary  mirror  before 
it  reaches  that  focus  where  a  knife  edge  was  located  in  the  conventional  sys¬ 
tem. 

The  wavefront  of  interest  can  then  be  reconstructed  by  illuminating  the 


hologram  with  the  reference  wavefront.  This  can  be  accomplished  by 
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developing  the  exposed  photographic  plate,  returning  it  to  its  original 

4 

position,  and  illuminating  it  with  only  the  original  reference  source. 

It  is  known  that  holographic  interferometry  can  provide  quantitative  in¬ 
formation  through  the  densi ty-foreign  relationship.  A  prerequisite  for  this 
application  however,  is  the  determination  of  whether  the  anticipated  refrac- 
tivity  change  falls  within  the  sensitivity  range  of  the  laser. 

The  refractive  index  of  a  fully  ionized  plasma  is  dominated  by  the 
electron  refractivi ty,  n.  The  fringe  displacement,  N,  due  to  the  substitu¬ 
tion  of  plasma  for  vacuum,  is  shown  below. 

L 

N  =  _/~(n-l )  dl/X  (1) 

o 

Here,  dl  is  the  element  of  path  length  along  a  ray,  L  is  the  distance 
across  the  test  area  cavity,  and  x  is  the  wavelength  of  the  laser  light  source. 
From  the  dispersion  relation  of  a  transverse  electromagnetic  wave  in  a  cold 
field  free  plasma  is 


K2  =  (w2  c2)/(l  -  wp2/w2)  (2) 

w  is  the  angular  frequency  of  the  laser  source  (w  =  2  t>  c/\)  and  w  is  the 

2  2 

cavity  gas  frequency  (w  =  4  v  a  c  /m).  Using  the  definition  of  the  index 
of  refraction  n  =  ck/w,  the  fringe  shift  relations  can  be  obtained.  For  the 
ruby  laser  (  x  =  6943°A) 
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N  =  oL/3.2  x  10 


(3) 
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_3 

o  is  the  electron  density  in  cm  .  The  length  L  is  also  in  cm.  Normal  dis- 

5 

cernible  fringe  displacements  are  between  1/10  and  30  fringe  per  cm.  De¬ 
tails  are  discussed  in  the  Appendix. 

As  a  matter  of  recording,  the  complex  amplitude  of  an  electromagnetic 
wave-front  possesses  both  magnitude  A(x)  and  phase  <j> ( x ) .  each  of  which  is 
essential  to  the  structure  of  the  wavefront.  The  magnitude  can  be  recorded 
by  simply  photographing  the  wave-front;  however,  the  phase  is  invariably  lost 
in  such  a  process  since  photographic  emulsions  are  generally  sensitive  only 
to  the  magnitude  or  absolute  value  of  the  scattered  amplitude.  Holographic 
recording,  on  the  other  hand,  provides  a  means  for  a  retrievable  recording  of 
the  phase  as  well  as  the  magnitude  of  the  field.  One  method  of  superposing  a 
coherent  reference  beam  on  the  field  striking  the  photographic  plate  is  as 
shown  in  Fig.  (3)  wherein  a  plane  wave  illuminates  a  region  containing  a  scat¬ 
tering  object  and  a  plane  mirror,  or  triangular  prism.  The  object  diffracts 
the  incident  radiation  to  generate  a  field  with  magnitude  A(x)  and  phase  0  (x) 
at  the  plate,  while  the  prism  turns  the  incident  plane-wave  through  a  small 
angle  0  to  form  a  field  with  a  uniform  magnitude  A  and  a  linear  phase  varia¬ 
tion  ax.  a  is  given  by  the  relation:^ 

a  =  2  tt  0  (4) 

The  equation  above  is  valid  for  small  angles  0  only.  Thus,  the  total  ampli¬ 
tude  striking  the  plate  is; 

A  eiuX  +  A{x)  ei9x 


(5) 
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Hence  the  intensity,  i.e.,  the  quantity  to  which  the  emulsion  is  sensitive 
is  given  by: 


I (x)  =  A2  +  A(x)2  +  2AqA(x)  cos  {a  x  +  <■  (x)) 


(6) 


As  indicated,  the  phase  of  the  scattered  wave-front  has  not  been  lost  in  com¬ 
puting  the  intensity,  as  it  would  have  been  if  the  reference  beam  were  not 
present.  The  emulsion  records  some  power  of  the  intensity,  that  is,  the 
amplitude  transmittance  of  the  resulting  photographic  plate  which  is  given 
by: 


T  (x)  = 


A2  + 
o 


A(x)2  +  2  AqA(x)  cos  [a  x  + 


(7) 


where  r  =  slope  of  the  Hurter-Driffield  curve.  The  Hurter-Driffield  Curve, 
shown  in  Fig.  (4),  relates  the  optical  density  to  the  exposure  of  the  holo¬ 
gram  for  a  particular  photographic  emulsion.  The  time  of  exposure  of  the 
hologram  may  be  related  to  the  intensities  of  the  object  and  reference  beams 
by  the  expression: 


E 


In  +  I„  +  2 


IRI0COS 


(Ky) ]  t 


(8) 


where  1^  =  Ref.  Beam  Intensity 

Iq  =  Object  Beam  Intensity 
t  =  Exposure  Time 

K  =  Constant  =  f  (incident  angles  of  object  and 
Ref.  waves). 
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The  optical  density  is  the  common  logarithm  of  the  opacity  8,  where 
2 

0  =  1/T  .  T  is  the  light  intensity  transmission,  i.e.. 


D  -  log-jQ  (6)  =  log10  (1/T2) 


=  -2  log1Q  T 


(9) 


The  linear  portion  of  the  H-D  curve  is  often  referred  to  as  the  range  for  cor¬ 
rect  tone  reproduction.  The  slope  of  this  curve  in  the  linear  range  is  the 
r  of  the  developed  emulsion.  For  a  given  material  the  value  of  r  can  be  con¬ 
trolled  by  changing  the  type  of  developer  and  the  time  of  development.7 

From  the  theory  of  optical  diffraction,  it  may  be  shown  that  the  complex 
amplitude  of  che  electric-field  vector  at  a  point  in  the  image-plane  of  an 
image-forming  system  is  the  Fourier  transform  of  the  distribution  of  the  com¬ 
plex  amplitude  of  the  electric-field  within  the  image  forming  aperture.  In 
other  words,  in  order  to  mathematically  obtain  the  distribution  of  the  complex 
amplitude  of  the  electric  field  vector  at  the  image  plane,  one  has  to  carry 
out  a  Fourier  transformation  for  each  point  in  the  diffraction  pattern. 
Mathematica ily,  the  above  facts  may  be  stated  as: 


\  <*». «)  ■  1  **&*>■//*. 

aperture 


exp 


iKA(x,y)j  exp  |-  iK(ax  +  3y) 


dxdy 

(10) 


Ep  is  the  complex  amplitude  of  the  electric-field  vector  at  a  point  P  in  the 
image-plane  in  the  vicinity  of  the  center  of  a  quasi -spherical  wave-front  of 
radius  R.  a  and  6  are  the  direction  cosines  defining  the  position  of  P  as 
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seen  from  the  aperture  centered  on  the  axis  of  the  quasi -spherical  wave-front, 
x,  y  are  the  coordinates  defining  the  aperture  A  (x,  y)  the  aberration  of  the 
wave-front  from  sphericity  and  k  =  2  /k where  A  is  the  wave-length  of  light 
used.  For  example,  in  the  case  of  a  rectangular  aperture, 

A/ 2 

exp  (-iKax)d  (11) 

-A/2 


Ip  (K.)  -  f 
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CHAPTER  IV 

Cold  Flow  Test  Program 

A  sketch  of  the  holographic  system  and  gas-  dynamics  facility  used  in 
the  study  is  shown  in  Figure  20.  The  components  of  the  system  are  listed 
below: 

a)  CR-5  Argon  Laser  (CW)  (Tunable  freq.)  Energy  is  5W  it  5145°A. 

b)  Korad  KQ1000H  Solid  State  Ruby  Pulse  Laser  -  6943°A  Energy  300  NW/cm 

c)  Spatial  Filters  XYZ  -  XY  Traverse  Stages 

d)  Beam  Splitters  -  80-20%  -  AR  Coated 

e)  Prism  Reflectors 

f)  Helium  Neon  Alignment  Laser  (15  mw  ) 

g)  6"  Spherical  Concave  Mirrors  -  60"  fl 

h)  Elliptical  Diagonal  Mirror  &  1"  by  3/4"  mirrors 

i)  Biconvex  Lenses  -  2"  and  4"  dia.,  30"  f.l.  and  15  "  f.l. 

j)  Schlieren  Knife-edge  XYZ  Traverse  Mount 

k)  Film-plate  Holder  and  Polaroid  Film  Holder 

l)  Polaroid  Land  Camera. 

The  system  was  designed  such  that  either  the  CW  Argon  laser  or  the  pulsed 
Ruby  laser  could  be  used  as  a  source  without  altering  their  position.  Light 
from  the  lasers  are  split  by  the  beam  splitter  in  an  80%-20%  ratio  such  that 
the  reference  beam  is  the  more  intense  of  the  two.  The  reference  beam  is 
then  reflected  by  the  prism  to  the  other  prism  on  the  opposite  table.  The 
object  beam  is  reflected  through  prisms  to  the  spatial  filter  S-j  where  it 
is  expanded  to  the  6"  spherical  mirror  M-| .  Thereafter  it  is  reflected  as 
a  collimated  beam  through  the  test  section  and  on  to  spherical  mirror  M2 
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that  directs  it  to  the  photoplate.  The  reference  beam  passes  through  the 
spatial  filter  that  expands  it  on  to  the  elliptical  diagonal  mirror  E. 

It  then  travels  through  the  30"  f.l.  lens  to  the  photoplate  to  form  the 
hologram.  When  the  Ruby  laser  is  used  the  light  travels  through  an  additional 
prism  reflector  on  to  the  first  beam-splitter  and  then  undergoes  the  same 
process  as  described  in  the  case  of  the  Argon  laser.  A  shutter  is  used  to 
control  exposure-times  when  the  Argon  laser  is  being  used.  The  holographic 
system  is  set  on  two  large  tables  weighted  with  concrete  and  on  air  cushion¬ 
ing  provided  by  inflated  inner  tubes. 

The  selection  of  components  of  the  holographic  system  was  based  on  the 
type  of  configuration  that  was  to  be  incorporated.  Since  the  system  is  a 
transmission  holographic  system, the  number  of  components  is  larger  than  a 
comparable  reflection-type  system. 

The  CR-5  Tunable  Frequency  laser  was  selected  in  order  to  study  the 
light  absorption  characteristics  of  the  different  mixing-gases.  It  has  a 
coherence  length  of  about  10  m  and  can  be  tuned  to  lase  at  different  fre¬ 
quencies  in  the  blue-green  range  of  the  visible  spectrum.  The  power  output 
of  the  CR-5  is  about  5  watts  without  the  etalon  and  about  2  Watts  with  it 
installed.  The  Korad  K1000QH  Ruby  laser  is  of  the  pulse-type  and  has  a 
coherence  length  of  about  100  cm.  Its  pulse-width  is  of  the  order  of  a  10-30 
nanoseconds  and  operates  at  a  single  frequency  of  6943 *A  with  a  threshold 

p 

lasing  value  of  300  MW/cm  at  about  5  Kilovolts. 

The  45°  prisms  were  used  since  metal  coated  reflectors  would  be  damaged 

by  the  power  of  the  unexpanded  Ruby  laser  beam.  Similarly,  the  80-20%  beam 

splitters  for  the  Ruby  laser  beam  were  dielectrically  coated.  The  spatial 
filtering  arran agements  for  each  of  the  lasers  was  also  different.  In  the 
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case  of  the  Argon  laser  the  beam  is  first  incident  on  a  microscope  objective 
(10X  or  20X)  before  striking  the  pinhole.  In  the  case  of  the  Ruby  laser 
the  beam  is  first  converged  by  a  simple  converging  lens  and  then  brought 
to  focus  at  the  pinhole  from  where  it  is  expanded  on  to  mirrors.  An  elliptical 
diagonal  mirror  was  incorporated  at  45°  to  the  reference  beam  before  passing 
through  the  bi-convex  lens  L^,  in  order  to  maintain  a  circular  beam  cross- 
section  as  it  rotates  through  90°  and  onto  the  photoplate.  The  order  of 
flatness  of  all  components  was  about  a  quarter  wave-length  of  the  Argon 
green  frequency.  The  diameter  of  the  6"  mirrors  was  selected  in  order  to 
cover  completely  the  test  section  windows  that  were  4x4  inches.  Also 
their  focal  length  of  60  inches  was  selected  because  the  longer  the  focal 
length  the  better  the  sensitivity  and  accuracy. 

The  reconstruction  system  employed  a  simple  15  biconvex  lens  to  bring 
the  reconstructed  wavefront  to  a  focus.  A  Schlieren  knife-edge  holder  with 
XYZ  traverse  positioned  at  the  focus  of  the  15  lens  enables  Schlieren  pictures 
to  be  made  using  a  Polaroid  film-pack. 

In  addition  to  the  TEMqq  mode  requirement  for  holography,  there  are 
other  important  considerations  necessary  in  laser  performance.  The  various 
modes  of  operation  of  a  laser  are  generally  not  coherent  with  each  other. 

The  lowest  order  mode  is  most  widely  used  because  the  Gaussian  intensity 
distribution  is  less  divergent,  more  easily  focused  and  is  most  easily  used 
in  calculations.  Population  inversions  between  several  pairs  of  energy 
states  allow  most  lasers  to  produce  more  than  one  frequency  of  light.  Such 
a  phenomenon  is  known  as  "frequency  broadening",  and  is  a  hindering  factor 
in  coherent-light-imaging  processes  such  as  holography.  In  order  to  obviate 
this  problem  it  was  necessary  to  suppress  other  undesirable  frequencies. 
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A  technique  that  could  be  used  to  accomplish  this  involves  turning  the 

reflector  inside  the  laser  cavity.  Another  simple  technique  is  the  use  of 

the  tilted  Fabry-Perot  etalon  which  is  inserted  into  the  laser  cavity.  The 

4 

transmissivity  of  the  etalon  is  given  by  the  relation  , 


2 

T  =  [1  +  4--Slp-  (2irnh  cosG'/X')]*1 

P  1  -K 


(12) 


This  reduces  cavity  efficiency  at  frequencies  other  than  the  etalon  resonant 
frequencies.  Figure  4  illustrates  the  procedure  just  described.  The  trans¬ 
missivity  maxima  of  the  etalon  are  matched  with  one  of  the  lasing  frequencies  by 
varying  the  angle  0’  such  that  the  peak  transmission  frequencies  are  shifted. 

The  approximation  is  good  when  6*  =  0,  which  was  generally  the  case.  Normally, 
the  laser  is  operated  at  the  second  peak  power  shown  in  Figure  4  to  realize 
full  suppression  of  undesired  modes. 

For  satisfactory  performance  of  the  optics  system  the  first  criterion 
was  a  vibration-free  mounting  of  components.  This  was  accomplished  by  fasten¬ 
ing  the  components  to  basemounts  with  rigid  3/4"  brass  rods  at  the  minimum 
height  necessary.  The  axis  of  the  Ruby  laser  was  6.688"  from  the  table 
surface.  All  optical  components  such  as  mirrors,  lenses,  beam-splitters, 
etc.  were  mounted  with  their  axis  of  symmetry  at  this  height.  The  tables 
were  mounted  with  their  tops  resting  on  air  cushions  provided  by  inflated 
inner  tubes.  The  pinholes  in  the  spatial  filters  were  critical  to  provide 
proper  alignment  and  diffraction.  Intensity  of  the  expanded  beam  was  found 
to  be  nonuniform  on  occasions  and  in  some  cases  insufficient  at  the  photo- 
plate.  This  was  remedied  by  using  pinholes  of  18  to  23  microns  in  the  object 
and  reference  beams.  It  was  found  that  the  quality  of  the  beam  greatly 
improved  with  an  increase  in  the  pinhole  diameter. 
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The  matching  of  path  lengths  prior  to  running  of  tests  was  carried  out 
by  individually  fixing  the  distances  between  various  components.  Certain 
values  were  automatically  fixed  by  optical  constraints,  e.g.  the  distance 
between  the  pinhole  from  spatial-filter  S-j  to  the  collimating  spherical 
mirror  was  fixed  by  the  focal  length  of  the  mirror  since  the  pinhole  is 
required  to  be  located  at  the  focus  for  perfect  collimation  after  reflection. 
Once  the  position  of  and  were  fixed  (assuming  that  is  centered  along 
the  axis  of  the  test-section)  the  angle  between  the  reflected  object  beam  and 
its  direction  before  reflection  (<j>j)  was  also  fixed.  Thereafter  the  angle 
<J>2  between  the  object  bean  directions  prior  to  and  after  reflection  of 
spherical  mirror  M,,  was  made  equal  to  in  order  to  ensure  that  the  object 
beam  was  incident  on  the  photoplate  parallel  to  its  direction  prior  to 
striking  collimating  mirror  .  The  orientation  of  the  photoplate  was  also 
adjusted  so  that  it  normally  bisected  the  angles  between  the  object  and 
reference  beams  at  the  plate.  These  details  are  indicated  in  Figure  18.  By 
positioning  the  lasers  and  reflectors,  the  paths  of  the  reference  and  object 
beam  were  adjusted  to  be  equal. 

To  properly  expose  the  plate  an  illumination  of  10  ergs/cm  for  a 
shutter  exposure  of  1  second  for  the  AGFA  10  E56  film  was  required.  A  cor¬ 
relation  was  obtained  for  photoplate-power  values  at  various  values  of  laser- 
current.  On  measuring  the  losses  in  intensity  through  various  components  of 
the  system  it  was  noted  that  the  prisms  transmitted  almost  99%  of  the  light 
while  the  mirrors  transferred  nearly  100%  of  the  incident  intensity.  However, 
it  was  observed  that  the  distance  between  the  spatial  filter  S,,  and  the  ellip¬ 
tical  diagonal  mirror  was  significant  in  that  it  determined  the  amount  of 
the  light  intensity  in  the  reference  beam  that  actually  reached  the  photoplate 


23 


via  the  elliptical  diagonal  mirror.  The  distance  between  the  two  was  adjusted 
such  that  a  maximum  possible  portion  of  the  expanded  reference  beam  would 
be  incident  on  the  elliptic  diagonal  mirror  and  hence  onto  the  photoplate. 

The  same  occurred  with  regard  to  the  distances  between  the  elliptical  diagonal 
mirror  and  lens  Lp  and  and  the  photoplate.  An  alternate  arrangement 
of  switching  the  position  of  the  spatial  filter  and  the  photoplate  was 
also  tried  successfully  to  improve  the  intensity.  These  alterations  are 
indicated  in  Figures  22,  23,  24,  and  25. 

The  alignment  of  the  Argon  and  Ruby  lasers  were  a  critical  and  time- 
consuming  task.  The  height  of  the  Ruby  laser  was  selected  as  the  standard 
beam  height  and  the  Argon  laser  was  elevated  to  that  height.  Alignment  of 
the  Argon  laser  was  performed  by  adjusting  the  heights  of  various  components 
in  the  optics.  The  CW  Argon  beam  was  convenient  to  use  to  align  that  system, 
while  for  the  Ruby  laser  a  CW  helium-neon  laser  was  positioned  alongside 
it  (shown  in  Figures  19,  20)  where  the  beam  was  directed  along  the  axis  of 
the  Ruby  laser  until  it  was  incident  on  the  beam  splitter  and  thence  to 
the  optical  set  up  for  the  Ruby  laser.  An  autocollimator  was  also  used 
to  line  up  the  individual  components  of  the  Ruby  laser  such  as  the  laser 
head,  Q-switch,  etalon,  rear  reflector,  etc.)  prior  to  lining  it  up  with 
the  optics  proper. 

In  order  to  test  the  stability  of  the  system  a  Michelson's  interfero¬ 
meter  configuration  (Fig.  21)  was  set  up  utilizing  the  Argon  laser  and  the 
beam-splitters  of  the  current  system.  The  stability  of  the  system  was 
checked  by  the  movement  of  the  interference  fringes  when  the  system  was 
slightly  disturbed.  It  was  noted  that  the  system  was  extremely  stable 
since  the  number  of  fringes  moving  across  the  field  was  negligible.  However, 
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the  system  was  not  entirely  insensitive  to  movement  as  the  table  tops  were 
not  rigid  but  were  freely  supported  by  the  inflated  tubes  underneath  them. 
Effects  of  a  hot  air  draft  from  the  laser  power  supply  was  also  indicative 
of  the  sensitivity.  The  power  supply  had  to  be  reoriented  with  its  air- 
draft  facing  away  from  the  system.  On  observing  the  object  beam  on  the  exit- 
side  of  the  test  section  it  was  observed  that  oblique  fringes  were  present 
across  the  square  cross-section.  These  may  have  been  caused  by  interference 
due  to  slightly  non-parallel  test-section  windows,  and/or  interference  due 
to  asymmetry  between  the  test-section  axis  and  beam-axis.  Adjustment  of 
the  test-section  windows  produced  a  slight  improvement  in  the  beam  quality 
on  the  exit-side  of  the  test-section. 

Based  on  readings  of  photoplate  illumination  at  various  laser-current 

o 

vclues  obtained  from  the  experiment,  (Figure  6)  and  a  value  of  10  ergs/cm 
required  for  an  exposure,  the  shutter  speed  for  the  Argon  laser  for  the 
various  tests  were  determined.  The  beam  diameter  in  each  case  was  assumed 
to  be  4".  The  times  for  developing  and  fixing  the  films  were  about  5  minutes 
each.  Developing  and  fixing  times  were  varied  for  double-exposure  tests, 
depending  on  the  extents  of  exposure.  The  final  image  is  obtained  by  re¬ 
constructing  the  flow  field  from  the  hologram  where  the  quality  and  efficiency 
of  reconstruction  is  extremely  important! 

The  efficiency  of  reconstructi on  n  is  given  by 


n  = 


Reconstructed  wavefront  intensity 
Reference  beam  intensity 


(13) 


where 


fl  = 


4R 

(l+R)2 


(14) 
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R  _  t  / t  =  Reference  beam  intensity 
1R/1o  ”  Object  beam  intensity 

£0  •  dc  exposure  term  ■  (I„  +  I0)tejq>  tjn]e 

6  =  slope  of  transfer  characteri Stic  curve  (plot  of  trans¬ 

mission  of  hologram  vs.  exposure) 

It  may  be  seen  that  for  a  fixed  E 

0 

_ R _ 

n  00  (HR)2 

which  indicates  that  for  a  low  value  of  R,  the  value  of  the  reconstruction 
efficiency  is  high.  However  a  low  value  of  R  also  increases  the  probability 
of  strong  higher  order  diffractions  that  may  severely  degrade  the  recon¬ 
structed  image.  As  a  compromise  between  efficiency  and  image  quality  the 
value  of  R  was  chosen  to  be  approximately  3  for  most  tests.  Experimentally, 
it  was  found  that  the  angle  between  the  object  and  reference  beams  signifi¬ 
cantly  influenced  the  resolution  of  the  holograms.  An  angle  of  10  degrees 
was  finally  chosen  for  the  best  resolution.  Also,  as  the  angle  between  the 
reference  and  object  beams  were  increased,  the  chances  of  spurious  reflec¬ 
tions  was  reduced.  It  was  noted  that  the  image  quality  of  the  hologram  was 
very  sensitive  to  its  alignment  with  respect  to  the  reference  beam.  A  slight 
rotation  about  a  horizontal  or  vertical  axis  was  found  to  cause  noticeable 
changes  in  the  clarity  of  images  which  were  near  the  diffraction  limits. 
Shutter  exposure  for  the  reconstruction  using  the  Polaroid  film-pack  was 
found  to  be  optimum  at  about  3-10  msec. 

Shown  in  Table  I  is  a  compilation  of  results  obtained  from  various  tests. 
The  aim  of  these  tests  was  to  establish  the  credibility  of  the  system.  Three 
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types  of  tests  were  performed: 

a)  Object-Reflection  Holograms 

b)  Gas  Flow  Tests 

c)  Flame- induced  density  gradient  tests 

Object  Reflection  Holograms  The  procedure  followed  was  to  remove  the  mirror 
and  replace  it  with  the  desired  object  for  record.  In  most  cases  both 
the  object  and  reference  beams  were  used.  In  one  case  only  the  object  beam 
was  used  to  form  the  hologram  by  placing  the  object  in  front  of  the  mirror 
and  blocking  off  the  reference  beam  while  exposing  the  plate.  In  another 
instance  the  reference  and  object  beams  were  switched  and  distances  between 
objects  and  photoplate  were  varied. 

Flame-Induced  Density  Gradient  Tests  These  tests  were  made  with  a  flame 
that  was  held  up  in  front  of  the  test  section.  Using  a  double  exposure, 
holograms  were  made  and  the  reconstructions  obtained  provided  high  quality 
interferograms  of  the  flow  around  the  flame. 

Gas-Flow  Tests  The  primary  gas  flow  tests  used  Nitrogen  and  Helium  as  the 
two  non-reacting  mixing  consti tuents .  Tests  showed  poor  quality  reconstruc¬ 
tions  partly  caused  by  a  temporary  laser  malfunction  but  also  because  of 
an  extensive  amount  of  spumous  reflections  from  prisms,  the  test  section, 
and  beam  splitters.  Blocking  off  these  increased  the  quality  of  the  holograms. 
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TABLE  II 


PHOTOPLATE  POWER-METER  READINGS 
(WITH  PLEXIGLASS  TEST  SECTION  WINDOWS) 


S.  No. 

Laser  Current 
(Amps) 

Ref.  Beam 
(watts) 

0^3 
T3  ft) 
a>  (jo 

13  05  CO 

CO  3  C 

0>  CO  -S 

CO  — <•  — <• 

O  *-H  3 
3T  3  fl> 

Cl>  CO  QJ 

1. 

12 

.01 

2. 

15 

.01 

CQ  3  C 

o>  c/>  -s 

CO 

3. 

18 

.012 

-j.  -..to 

-*•  -i.CQ 

4. 

21 

.014 

O  c+  c+ 

£  O  0> 

0)  -5 

0>  3 
~o  fD 

O  rt  r+ 

^  o  0> 
0)  -s 
-1 

5. 

24 

.015 

6. 

27 

.02 

7. 

28 

8. 

30 

0.005 

TOTAL 

Ref.  +  Obj .  Beam 
(watts) 


.012 

I 

.012 

'  ■ 

.010 

1 

.012 

£H 

.014 

IS 

.014 

.021 

.019 

.016 

.021 

.020 

Inferences:  Power  meter  either  insensitive  to  small  power  changes  or 
was  malfunctionary. 
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Interference  Grating 
(  Carrier) 


Amplitude  and  Phase  Modulated 
Interference  Grating. 


FIG.  3  .  Amplitude  and  Phase-Modulated  Interference-Grating. 


Photoplate 


FIG. 4  Recording  of  hologram  for  a 
plane-wave  generating  object  (Ml) 


FIG.  5.  Recording  of  halagram  of  3-d  object  (  Stroke  G.W. 


FIG.  Interference  pattern  of  a  diffraction  grating 
7.  illuminated  by  a  plane-wave 
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FIG,  8,  Spherical  Wavefront  Interference  pattern 
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FIG.  9 
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.  Fourier  Transform  Reconstruction  from  Fourier-  Transfc 
Hologram  H(x).  (  Ref.  6  ) 
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FIG.  10.  Hologram-  Recording  Parameters. (  Ref.  6  ) 
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FIG.  11 .  A  Typical  Hurter-Drif field  Curve  for 

a  photographic  recording  medium, 

(  Reft  7  ) 


FIG,  12.  Variation  of  Reconstruction  efficiency^  for 

fixed  p  and  E  values  as  a  function  of  R=  I_/I. 

'  O  R  0 

(  Ref.  7  ) 
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FIG. I  3  .  Typical  Transfer  Characteristic 
Curve  for  a  photographic  medium 
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14 ,  A  simple  method  of  recording  a  hologram  by 
superposition  of  a  coherent  background  over 
a  scattered  field.  (  Ref.  4  ) 
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FIG.  15.  The  laser  Cavity  (  Ref  4  ). 
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FIG.  16*  Possible  Resonant  Frequencies  Centered  about  one 
Colour  (  Ref.  4  ). 
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[FIG .  1 7~  The  tilted  etalon  in  a  Laser  Cavity. 

(  Ref.  4  ). 
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Etalon  at  6  +  6' 


Resonant  Frequencies 
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FIG. _[8,  Tilted  Etalon  Characteristic  Curve.  (  Ref.  4 


Note  i"At  non-resonant  frequencies  of  etalon  losses  are  introduc 
Lasing  action  becomes  frequency  dependent.  Proper  selecti 
and  adjustment  of  etalon  (e.g.  tilting)  sharpens  the 
resonance  curve,  thus  tuning  laser  to  single-frequency 
mode.  Tilting  an  etalon  by  more  than  a  fraction  of  a  mill 


radian  causes  severe  losses  even  at  resonance."  (  Ref.  4 
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FIG,  19 .  Spatial  Filtering  Arrangement  for 


Lasers.  (  Ref  (1)  ). 
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Intensity  Distbn,  at 
A  for  Gaussian  Input 


Note  i  In  the  case  of  the  Argon  laser  the  simple  lens  L  is 


replaced  by  a  compound  -lens  system,  usually  a  micro¬ 
scope  objective  of  20  X  or  10  X. 


THROTTLH  VALVE 


FIGURE  20.  HOLOGRAPHIC  AN'D  GAS  DYNAMIC  FACILITY 


qure  22. 


h-18.46—1 


gure23.  Nozzle  Design  No.  1  (dimensions  in 


0-032  (VHtoAf) 


No.  2  (from  "LAMP"  Program). 


Table  III.  Coordinates  of  Nozzle  No. 


CO-ORD ■  SYSTEM  hOR  NOZZLES  w ITH  0.0  ,\i  C£FR£  Oh  Each 


NOZZLE  -  Cl  V  SCALED  1:6-55%  ASSEMBLY  all  in  in a-t t  S 
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ANALYSIS 

Tests  indicated  that  the  contoured  nozzles  1,  2  and  4  (Figures  22,  23, 

24,  25,  26,  and  30)  provided  better  control  and  repeatabi 1 i ty  than  the 
straight  nozzle  3  (Figure  27  and  28),  although  the  straight  nozzle  provided 
occasional  high  mixing  levels.  It  is  expected  that  the  favorable  pressure 
gradient  created  through  the  contour  and  the  stabilized  boundary  layer  on 
the  internal  walls  of  the  contoured  nozzles  provided  the  consistent  flow 
conditions  downstream.  The  less  stable  boundary  layer  internal  to  the 
straight  nozzle  increased  the  chance  of  separation,  adverse  compression 
shocks  and  thus  variable  and  unpredictable  flow  conditions  downstream.  The 
alternate  large  and  small  nozzle  arrangement  was  not  found  to  greatly  influ¬ 
ence  mixing.  Data  and  pictures  suggested  that  the  comparison  of  flow  field 
mixing  downstream  of  nozzles  1,  2  and  4  were  not  particularly  sensitive  to 
geometry. 

The  study  concerning  density  mapping  was  only  partially  successful  in 
that  the  procedure  was  identified  and  applied  to  general  flow  fields;  however, 

applications  to  conditions  of  this  study  was  hampered  by  the  loss  of  the 
Argon  laser  at  a  critical  point  in  the  study.  Quantitative  data  was  there¬ 
fore  not  obtained  as  to  the  evaluation  and  diagnosis  of  flow  field  mixing. 

Qualitatively,  test  nozzle  4  provided  much  information  on  mixing.  For 
example,  the  high  mass  flow  of  an^  He  mixture  in  the  center  nozzle  (Figure 
31),  which  provided  a  lower  molecular  weight  and  higher  Mach  number  there, 
reduced  the  mixing  significantly  in  the  near  field.  Density  level  variations 
of  the  same  flow  conditions  (observed  in  the  Figure  32  interferogram)  also 
indicates  low  level  mixing.  The  sensitivity  of  flow  conditions  is  observed 
in  Figure  33,  where  good  mixing  is  occurring.  The  only  significant  difference 
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in  the  flow  conditions  from  Figures  31  and  32  was  the  exit  pressure  conditions 
for  the  center  nozzle  (which  was  lowered)  that  forced  rapid  inward  flow  from 
the  outer  nozzles.  Figure  34  is  an  indication  of  moderate  to  poor  mixing. 

Here  the  mass  flow  of  adjacent  nozzles  were  nearly  the  same. 

As  a  base.  Figure  35  was  run  with  near  zero  flow  rate  through  the 
center  nozzle  and  methane  CH^  on  either  side  at  21  grams/sec,  which  should 
have  maximized  the  divergence  of  the  outer  nozzles  in  the  near  field,  yet 
little  of  this  occurred.  A  similar  flow  pattern  occurred  for  matching  all 
flow  properties  for  the  nozzles  (Figure  36).  Figures  37  and  38  are  hori¬ 
zontal  and  vertical  (knife  edge)  schlierens  of  high  center  mass  flow  rates, 
low  exit  pressures  of  the  center  nozzle  and  a  relatively  high  cavity  pres¬ 
sure,  providing  moderate  to  good  mixing.  Lowering  the  cavity  pressure, 
matching  the  flow  rates,  but  increasing  the  exit  pressure  differences  pro¬ 
vided  very  good  diffusion  mixing  (Figure  39),  which  also  reduced  the  shock 
strength  and  viscous  effects. 

As  a  check  on  the  effect  of  maximizing  the  flow  rate  of  the  exterior 
nozzles.  Figure  40  provides  additional  evidence  that  thi s condi tion  does  not 
enhance  mixing.  Conversely,  increasing  the  flow  rate  of  the  center  signifi¬ 
cantly,  plus  a  higher  cavity  pressure  (about  3  times)  provided  excellent 
mixing  of  the  flow  streams  from  the  nozzles  (Figure  40). 

Demonstrated  in  the  tests  was  the  mean  bulk  flow  direction  change  toward 
the  nozzle  that  provided  the  least  momentum  flux  and  was  typical  for  all  the 
runs  and  dominant  in  the  non-equilibrium  zone  of  the  near  field,  which  occurred 
up  to  10  diameters  downstream.  Solutions  for  the  non-equilibrium  zone,  assuming 
a  laminar  mixing  layer  was  obtained  from  the  momentum  integral  method  and  pro¬ 
vided  approximations  and  insight  to  the  mixing  zone,  layer  thickness  and  length. 
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with  reasonable  agreement  with  tests.  Approximations  were  also  obtained 
for  the  turbulent  mixing  layer  length,  suggesting  the  test  conditions  were 
mostly  transitional  to  turbulent  in  the  non-equilibrium  zone. 
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Figure  35.  Shadowgraph  of  Test  Flow  (Zero  Center  Flow  -  Run  9). 


Figure  36.  Shadowgraph  of  Test  Flow  (Matched  How  -  Run  11). 


Appendix  A 

Shadow,  Schlieren  and  Holographic  Fundamentals 
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The  character  of  gas  dynamics  makes  it  possible  to  view  turbu¬ 
lence,  shock  patterns,  boundary  layers  and  other  phenomena  in  the 
flow,  caused  by  variations  in  properties,  particularly  density,  in 
several  ways.  This  is  possible  because  of  the  manner  in  which  light 
passes  through  the  medium  being  investigated.  The  figures  below 
illustrate  how  light  may  be  refracted  through  a  gaseous  medium,  where 
the  density,  p,  increases  in  the  positive  y  direction. 

The  index  of  refraction  n  =  c0/c  =  1  +■  p«/:s  . 

The  angle  of  refraction  c  =  /*"  d<t» 

o 


From  geometry,  dn  =  dy  sin  a,  and  d?  = ; dz  sin  'i.  And  since  c  =  cQ/n, 

we  can  write  (dc)  -  c  dn 

dy  n  dy 

where  d4  =  tc ,  the  distance  traveled  along  r^. 


The  angle  of  refraction  can  thus  be  written: 


j  dr 


f  =  f _ 

J  d  J  dy  sin 

/C  dn\  =  f  dr- 
Ct  V  n  dy  /  J  sin 


I  l 


_L  dn 
n  dv 


1. 


L 


/ 


_1  dn 
n  dy 


dz 


dz 


From  a  three  dimensional  view  point,  as  shown  in  Tigures  below,  the 
direction  of  refracted  light  will  depend  on  the  direction 
flow.  For  the  case  below  the  flow  is  the  x  direction,  while  the 
light  beam  is  in  the  z  direction.  The  resulting  angles  of  refraction 
in  the  y  and  x  directions  can  thus  be  written. 


TEST  SECTION 


RAY  OF  LIGHT 


The  shadow  method  for  making  shadowgraphs  is  illustrated  below 
Here,  the  light  is  refracted  as  it  passes  through  the  test  section, 
providing  light  and  dark  regions  on  the  photo  plate  or  screen.  The 


shadow  effect  (dark  and  bright  regions )  depend  on  the  relationship. 


OR  SCREEN 


At  the  test  section  ^  =  Lh  df 

«s  dy 


The  illumination  I  =  3t/:iy 

We  may  therefore  write  in  plane  flow 
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Schlieren  photographs  are  characteri zed  by  symmetric  light  and 
dark  regions  caused  by  density  gradients  in  the  direction  of  the 
light.  The  figure  below  illustrates  a  Schlieren  system  for  pencils 
of  light  as  they  proceed  through  a  test  section  of  uniform  density. 

A  knife  edge  placed  at  the  focal  point  intercepts  the  light.  At  the 
knife  edge,  a  portion  of  the  light  is  refracted  by  a  distance  -\y=f  s 
The  contrast  illuminated  on  the  photo  plate  is 


Q  ~  I  1 


where,  c  --  •'■v /y,  -f  e  /y  ,  and  e  -  J  ^ 

y  ”,  ay 

we  can  thus  write  a  relationship  depicting  the  contrast 


THE  KNIFE  EDGE 
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The  fundamental  operations  of  the  interference  method,  are 
illustrated  in  the  figures  below.  As  shown,  light  travels  to  a  beam 
splitter  (BSI),  generally  passing  50%  of  the  light.  The  separated 
light  then  travels  to  respective  mirrors,  each  traveling  through 
different  media.  The  light  is  then  gathered  at  the  second  beam 
splitter  (BS2),  each  light  ray  possessing  a  different  character  be¬ 
cause  of  their  now  different  optical  paths.  The  light  finally  passes 


through  a  lens  and  onto  a  photo  plate  that  illuminates  the  fringe 
patterns  caused  by  the  different  wave  character.  The  fringe  and  den¬ 
sity  relation  for  the  interference  method  is  shown  below. 


ROTATION  OF  BS2 
AN  ANGLE  t 


6;-. 


Fringe  ami  iVnsi  tv  i<c:  In  t  ion 


N 
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Appendix  B 


Density-Refractive  Index  Relationship 
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Density-Refractive  index  Relationship 

There  are  three  impirical  expressions  that  have  been  proposed  for  the 
relationship  between  the  refractive  index  n.  for  a  particular  wavelength  of 
radiation,  and  the  density,  p,  of  a  given  transparent  body.  These  expressions, 
in  historical  order  are: 

2 

Newton-Laplace  (n  -l)/p  =  Constant 

Gladstone-Dale  (n-l)/p  =  Constant 

2  2 

Lorenz-Lorentz  (n  -l)/(n  +2)p  =  Constant 

The  first  expression  was  developed  by  Newton  and  gave  a  reasonable  value 
for  most  liquids  when  the  temperature  was  varied,  and  was  partially  based  on 
the  emission  theory  of  light.  Laplace  inferred  the  same  expression  from 
theoretical  considerations  based  on  the  form  of  the  undulatory  theory  accepted 
at  his  time.  As  a  result  of  experimentation  on  water  and  its  vapor,  it  was 
considered  that  the  relationship  did  not  hold  when  a  change  of  state  occurred. 

Better  agreement  with  experimental  facts  was  obtained  with  the  second  experi¬ 
mental  expression,  which  was  proposed  by  Gladstone  and  Dale  (1858,1863)  as 
a  result  of  experiments  with  water  and  carbon  disulphide  at  different  tempera¬ 
tures.  Sutterland  (1889)  attempted  to  give  a  theoretical  basis  for  this 
relation,  but  without  much  success. 

The  third  relation  was  independently  derived  by  Lorenz  (1880)  of 
Copenhagen  and  Lorentz  of  Amsterdam,  from  entirely  different  theoretical 
approaches.  Lorenz  considered  the  passage  of  light  through  media  consisting 
of  spherical  molecules  immersed  in  an  ether  having  the  properties  of  an 
elastic  solid.  Lorentz  derived  the  same  expression  from  his  own  theory  of 
electrons  considered  in  relation  to  Maxwells'  system  of  equations  for  electro¬ 
magnetic  theory  of  light.  Lorenz  and  others  have  found  it  particularly  suc¬ 
cessful  in  correlating  the  facts  when  a  liquid  passes  to  the  gaseous  state. 
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The  experimental  comparison  of  the  three  relations  in  the  case  of  gases 
is  not  so  easily  made  as  in  the  case  of  liquids,  although  a  considerable 
amount  of  research  has  been  devoted  to  the  subject  because  of  the  ease  with 
which  the  density  of  a  gas  may  be  varied  by  alterations  of  temperature  and 
pressure.  The  refractivity  of  gases  at  normal  densities  is  relatively  small. 

The  first  and  third  relations  then  reduce,  as  a  first  approximation  to  the 
Gladstone-Dale  equation. 

Interferogram  Redu ction 

The  final  result  of  any  of  the  six  techniques  of  holographic  interferometry 
is  an  interferogram  with  interference  fringes/ 

The  next  step  is  to  obtain  the  density  distribution  by  a  suitable  reduction 
from  the  interferogram.  Before  we  can  do  this,  we  have  to  establish  the  fringe- 
shift  density  relationship. 

The  following  are  a  few  basic  assumptions  made  from  geometrical  and  physi¬ 
cal  optics: 

1)  In  a  homogenous,  isotropic  medium  a  beam  of  light  may  be  analyzed  into 
rays  which  are  represented  either  as  vectors  or  straight  lines  in  3 
space . 

2)  A  monochromatic  ray  of  light  consists  of  a  sinusoidal  disturbance  of 
frequency  v  and  wavelength  X  propagating  along  the  normal  to  the  direc¬ 
tion  of  vibration  with  velocity  u. 

3)  X v  =  i)  is  satisfied  where  v  is  constant  and  A,u  are  functions  of  the 
coordinates . 

4)  Index  of  refraction  ’ n '  is  related  to  velocity  by  n  =  c/u  where  c  is 
velocity  of  light  in  vacuum. 

5)  For  gaseous  media  at  moderate  pressures  there  is  a  direct  relation 
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between  density  increments  p,  und  index  of  refraction  n-nQ,  thus, 

n-nrj  -  k(;  p.)  (1) 

where  k  is  the  G1 adstone-Dal e  constant. 

6)  For  axisyminelric  flues,  the  refraction  of  any  ray  in  the  disturbed 
beam  is  negligible. 

With  these  prel iminaries,  we  can  now  derive  the  fringe  shift  equation. 


v 


xo  dx  x 


Consider  a  general  density  disturbance  in  one  beam  of  a  Mach-Zehnder  inter¬ 
ferometer.  Coordinate  directions  are  defined  by  an  orthogonal  right-hand  triad 
with  x  parallel  to  the  interferometer  beam  and  y  vertical  as  in  the  Figure 


above. 


If  the  i-dex  of  refraction  function  n(x,y,z)  has  only  a  finite  number  of 

discontinuities,  the  fringe  shift  6(y,z)  is  given  by 

/<•  i  ,  w*  1 

l(y,z)  =  M  (n-n0  )dx  =  -*■(  (.-pjdx  (c 


)d.x  =  '  (.  -  pjdx 


Where  n  and  .  are  the  reference  index  of  refraction  and  density  in  the 

o  o 

uitisturbed  beam  and  >  is  the  vacuum  wavelength  of  light. 

"n  rove  (2),  we  observe  that  t.he  optical  pathlenqth  of  the  disturbed  ray 
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The  difference  Nd-NQ  gives  the  number  of  multiples  of  2f  by  which  the  two 
rays  differ  in  phase  when  recombined.  This  number  is  equal  to  the  fringe 
shift  f(y,z) . 


i  .e.  N 

o 


A1 

f  dx 

J  >.  ' 

X,  0 

o 


A1 

J 

x0 


n  dx 

_o_ 


Two-Dimensional  Density  Field s 


n-n 


If  n  =  n(y,z)  and  (■  -  e.)  =  -r — ,  then  for  two  dimensional  flues. 


*(y,z) 


- !  I’t, 


II  \A  2:ik 
;  u-pjdx.  v*- 


(3) 


2__k 

X* 


u-e.)[x] 


(4) 


the  flue  'is  in  the  y-z  plane,  the  term  ( ,  o )  can  be  taken  out  of  the  integral 
sign.  This  means  that,  the  density  distribution  is  two-dimensional  in  the 
y-z  plane. 


or 


c(u ,z) 
( o  -  (>* ) 


v^/(r-a)(x,-x0) 

fly.zh* 

2r7kt)q-x0T 


or 


o(y,z) 


=  (?• 


,  5(y,z)A» 

2nkXx]^xoy 


(5) 

(6) 


Thus  knowing  ( y , z )  from  the  i nterferogram  and  knowing  we  can  deter¬ 

mine  n(y,z) . 

Determ ination  of  F r i nge - S hj f ts  f rom  I  n_t  e rferograms 

To  obtain  the  graphical  plot  of  the  fringe  number  functions,  the  negative 
is  used  in  a  photo  enlarger  to  image  directly  upon  graph  paper.  With  the 
negative  adjusted  to  match  the  scale  of  the  paper,  the  positions  of  the  maximum 
and  minimum  of  the  figures  along  the  desired  cross-section  arc  determined 


visually.  Initial  fringe  numbering  is  arbitrary.  Commencing  with  a  fringe 
well  to  one  side  of  the  field  and  proceeding  across  the  field,  the  length 
requires  representing  integer  values  of  the  fringe  number. 


If  the  fringe  number  N  is  plotted  against  position  y,  a,  ab  issue,  a 
straight  line  is  obtained  for  fringes  in  the  undisturbed  field.  With  flow  in 
the  test  section  established,  the  plot  deviates  from  a  straight  line  in  the 
disturbed  part  of  the  field,  so  that  at  each  y  value  the  fringe  shift  is 
given  directly  by  the  difference  in  the  ordinalis  of  the  two  curves.  The 
fringe  shifts  are  actually  not  plotted  out  but  calculated  from  the  equivalent 
equation; 


My)  = 


I  N,-N 

NT|y(  d  0 


N  -N 

)  +  VYn(;rV-°) 


L  y1y0  1  °  °Vo 


(7) 


This  equation  is  obtained  as  follows: 


The  fringe  shift  <s(y)  at  any  point  a  distance  y  from  the  origin  is  given 

by 

f,{y)  =  N'-N 

where  N'  --  fringe  number  in  the  disturbed  region 
N  --  fringe  number  in  the  undisturbed  flow 


A 
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Now,  we  know  that  in  the  undisturbed  region  the  fringes  are  parallel  and 
equally  spaced.  Hence,  we  can  assume  N  to  be  linearly  related  to  y(  that  is 

N  =  Ay+B  (8) 

N  =  NQ@y  =  yQ 
N  =  N-j@y  =  y, 

Where  N-j  and  N^  are  two  fringes  close  together  in  the  undisturbed  region, 
close  to  the  boundary  of  disturbance 

y  and  y^  are  the  correct  positions  of  and  Nq  in  the  particular  cross- 
section  being  considered. 

M0  =  AVB  A  =  ,  B  =  V0(^> 

N-j  =  Ay-j+B  yly0 

5(y)  =  N  (Ay+B)  (9) 

The  sign  is  to  be  chosen  so  that  positive  fringe  shifts  result  from  retarded 
rays,  where  the  average  density  through  the  flow  is  greater  than  in  the  undis¬ 
turbed  region.  The  form  of  (9)  may  be  varied  to  suit  the  computing  machine. 

After  calculation  of  <$(y)  for  different  values  of  y  from  (9),  a  smooth 
curve  can  be  drawn  through  the  points  [y,<$(y)]. 


76 


Appendix  C 


Computer  Programs 
Including  Subroutines 


$ 


REAL  MU< 161),  MUPOS, MA,  MB*  MU  120* MU220,  M>US,  ME(  1 6 1 > 

DIMENSION  U<  161, 162>,ZU0<  161>,A<  162>,U1< 161) ,U2< 161),VES< 161) 
,  PD<  1 6  1  )  ,  P<  16 1 ) , DI  SP<  1  6  I  )  ,  AI  (  161),  TY  (  161),  TY  1  (  I  6 1 )  ,  T  (  1  6 1  > 

READ  <2#  15, END-9999)  (AI(OI),OI-l,  161,4) 

READ  <2,  1  I, END-9999)  MA, MB, CPA, CPB, CVA, CVB, AM, BM, TO, PO, 

MU220, MU  120, D2, D 1 
FORMAT  C6F10.4) 

FORMAT  (6E20.5) 

FORMAT  < IX, 161F10.0) 

FORMAT  (10F5.4) 

CALL  HEATR  (  M A,  MB,  CPA,  CPB,  CVA,  CVB,  AM,  BM,  R,  CM,  CP) 

CALL  THROAT  < R, CM, TO, PO, TS, PS, ZUOS, US ) 

CALL  VISCOS  ( TS, MB, MA, HU 120,MU220, Dt,D2, BM, AM, MUS ) 

VELO-US* *3048 
TK-<TS-460-32)*5./9*273 
ZU0G-ZU0S*453.59/<  30. 48)** 3. 

MUPOS-MUS* 14.882 
PTOR-PS*760/ 14.7 

WRITE  (3,12)  MA, MB, CM, R, TO, CP, VELO, TK, ZUOG , MUPOS, PTOR 

DX -.87  16/<  160.*. 1) 

D-. 1/12. 

RS-. I/2./I2* 

TSI -0.5 
DSI-TSI / 160. 

PD<  1)-PS*  144*32. 2/ CZ U05*US**2.  ) 

RES-ZUOS*  US*  D/MUS 
DO  1  J  1-2,  162 
U<  1,J1)»I. 

DO  2  11-1,161 
U< I  1,  l)-0. 

DISPCI 1 ) -0 • 

DO  16  01 1-3, 159,4 

A1(JI 1)»(A1<JI  l-2)*Al<0I  1*2) )/2. 

DO  17  OI2-2, 160, 2 

AI <0I2)-< Alt 012- 1>*AI <012* 1) )/2. 


DO  9  11-1,10 
NU-  I 

CALL  PRESS <  PS, R, PO, AI , DI SP, P) 
DO  6  0  11*1,17  1 
ZU0(0  1  !)•  1 . 

MU(0  1  D-MUS 
T<0  1  D-TS 


on  o  o  o  o 
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THIS  SUBPROGRAMME  IS  USED  TO  PREDICT  PROPERTIES 
AT  THROAT  FROM  SOME  RESERVOIR  CONDITIONS 

SUBROUTINE  THROAT  (GAMA*CM*TO*PO*TS*PS*ZUOS*US) 

TO-TO*9/5*32*460 
ZUOO-PO* 1 44* < CM/ 1544) /TO 
TS-2./<GAMA+ 1 .)*TO 

PS-( 2. /(GAMA* !•)>**< GAMA/C  GAMA- 1 ) )*PO 
ZUOS-< PS/PO>**< 1/GAMA)*ZUQ0 
US- SORT ( GAMA* 1 544 /CM*  32 . 2*TS> 

RETURN 

END 


THIS  SUBPROGRAMME  IS  USED  TO  DETERMINE 

SPECIFIC  HEAT  RATIO  AND  MOLECULAR  WEIGHT  OF  MIXED  GASES 
SUBROUTINE  HEATR  ( MA*MB*  CPA*  CPB*  CVA*  CVB*  AM/ BM* R*  CM#  CP> 
REAL  MA*MB 

CP-<MA*CPA*MB*CP8)/<MA*MB> 

CV-  (  MA*CVA+MB* CVB>  /<  MA*MB> 

R-CP/CV 

CM-C  MA*MB)  /(  MA/AM+MB/BM) 

RETURN 

END 


THIS  SUBPROGRAMME  IS  USED  TO  DETERMINE  VISCOSITY  OF 
MIXED  GASES  AT  TEMPERATURE  T 

SUBROUTINE  VISCOS  ( T* M2, M 1 , MU  1 20, MU220* D 1, D2* W2* V 1 * MU) 

REAL  Ml*  M2*  MU  1  *  MU 2*  MU 120*MU220*  MU 
T20C-528. 

DELMU1-<T-T20C>*D 1*5/9. 

DELMU2-<T-T20C>*D2*5/9. 

MU 1*MU 120+DELMU 1 
MU2-MU220*  DELMU2 
X1-(M1/V1)/<M1/V1*M2/V2> 

X2-<  M2/V2) /(M1/W1  +  M2/V2) 

IF  (X2.EQ.0.)  GO  TO  1 

PH  12-(  1*<MU1/MU2>**0.5*< V2/V1>**0. 25)* *2. 

PH  1 2- PH  12/(8  .**0.5)/<  1*W1/V2>**0.5 
PH2  l-< l*CMU2/MU 1)**0.5*C W 1/W2)**0.25)**2. 
PH21-PH2l/<8.**0.5)/<  1*W2/W  i>**0.5 
MU-MU  l/<  1  +  X2*PH  12/Xl)+H0fcM  U*l* PH21/X2) 

GO  TO  2 

1  MU-MU  1 

2  MU-MU*  6.7 1 9  7  E- 08 
RETURN 

END 
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DO  5  I-l.  160 
PR»P(  I*  1) 

PD  ( I  ♦  1  >  -  P(  I  ♦ 1>* 144.*32.2/(ZU0S*US**2. ) 

DO  3  J2- 1, 161 

3  A(J2)«ZU0(  J2)*U(I.J2)*MU(J2)/MUS 
DO  4  J«2.  161 

TERP-( - l./(ZUO(J)*U(I,J) >)*(PD(I*  1 )  -  PD  (  I  )  )/DX 
AA 1* ( A(J)+A(J+l))/2. 

AA2-(A( J)+A( J  -  1) )/2. 

TERMU«AAl*(U( I.J* 1) -U< I.J) ) -AA2*(U(  1  .  J) -U<  1  . J- 1) ) 
TERMU-TERMU/(DSI**2.*RES> 

4  U(I* l.J)«U(I.J)*DX*(TERP*TERMU> 

U(I  +  1.  162) *U(  I  ♦  1,  160) 

DO  8  J3«l. 161 
UI(J3)»U(I* 1.J3) 

T(J3)-TO-(Ul( J3)*US)**2./(2.*32.2*778.*CP) 
ZU0(J3)«P( l* 1>* 144. /( < 1544./CM)*T(J3) )/ZUOS 
MU( J3) «MUS*  (T(J3)/TS)**/75 
VES(J3)-SORT(R*  1  544/CM*  32. 2*T( J  3) > /US 
8  CONTINUE 

CALL  DISPLE  ( D. ZUO. DSI . U 1 . VES. TY. TTY, TY 1 , ME. DI SP2) 
IF  ( I /EQ* 160)  GO  TO  50 
DISPU+  D-DISP2 
50  CONTINUE 

IF  (II.LT.3)  GO  TO  5 
IF  <1. NE.NW)  GO  TO  5 
NW-NW+40 
DXl-DX*. 1/12. 

DX2-DX 1*30.48 
RS 1-RS*  30. 48 
X 1  *  DX  1*1 
X2»X 1*30.48 
XDR«X 1/RS 

PTOR«P( I*l)*(76 0/14.7) 

DPDX-(  P( I* 1) -P( I ) ) / <  DX  * • 1/ 12)* <  760/ 14.7) /30. 48 
WRITE  (3. 14)  II. XI.  X2.  RSI. TTY 
WRITE  (  3#  20)  I  .XDR.  PTOR.  DPDX.  DX  1*  DX2 
DO  30  K-1.20 

YDR-AI ( I  ♦ l)/.05-TY(K)/(R5*  12.) 

VELO-U 1(K)*US* .3048 

TK»(T(K>-460-32)*5./9*273 

ZUOG-ZUO(K)*ZUOS*453.59/<30.48)**3. 

MUP05-MU <  K ) *  14.882 

30  WRITE  (3.21)  K. YDR. VELO. TK.Z UOG. ME( K) . MUPOS 

5  CONTINUE 

ERROR- ABS< (DISP< 16 1 ) -DI SP2) /DI SP2) 

IF  (ERROR.LT. 0.001)  GO  TO  9999 
DI  SP(  16D-DISP2 


9  CONTINUE 


14 

FORMAT 

(/1X.I4.4E20.7) 

20 

FORMAT 

(/IX. 14. 5  E2  0.7) 

21 

FORMAT 

(  1X.I3.6E18.5) 

9999 

STOP 

END 

THIS  SUBPROGRAM  IS  USED  TO  DETERMINE  POESSURE  DISTRICUTION 
ALONG  THE  NOZZLE  AXIS 

SUBROUTINE  PRESS  CPS,R,PQ,A,  DISP,P> 

DIMENSION  PC  161), AC 161),DISPC 161) 

PC  1>»PS 
DO  I  1-1,160 

ARAT-C  ACI+  D-DISPC  1+  l))/.05 
CALL  EX  I T  PC  R,  PO, ARAT,  PX ) 

PCI*  1)-PX 

RETURN 

END 


THIS  SUBPROGRAMME  IS  USED  TO  DETERMINE  DISPLACEMENT  THICKNESS 
SUBROUTINE  DISPLE  C  D, ZUO,  DSI  E,  U,  VES,  TY,  TTY, TY  1  ,ME,  DI  SP) 

REAL  MEC  17  1) 

DIMENSION  UC 161>,DELYC 161),DELY1C 161), TYC 161>, 

TY 1 C  16 1 > , ZUOC 161), VESC 161) 

TYC l>-0. 

DISP-0. 

TY  1  C  160)-. 2936 
DO  8  J  1-160,2,-  l 

DELY  1C J  D-DSIE/CCZUOCJ  D+ZUOCJ  1+  1)  )/2*CUCJ  |)*UCJI+  1 )  ) /2  .  )  *  1 2  .*D 
TY1CJ1-D-TY1CJ1)  -DELY  1CJ  1) 

DO  7  J-l, 160 
MECJ)-UCJ)/VESCJ) 

DELYC J )-DSIE/CCZUQCJ)*ZU0CJ*l))/2*CUcK)*UCJ*  1)  )/2.)*  12*D 
I r  CUC J) .EO.UC J* 1) )  GO  TO  9 
U1-UCJ)/U<  161)*ZU0CJ)/ZU0C  160) 

U2-UCJ+D/UC  161)*ZU0CJ)/ZU0C  160) 

DISP-DISP*C<  l*Ul)t(  1-U2)  )  /2.*DELYCU) 

TTY-TYC J ) 

TYC  J* 1 ) -TY  C  J )  *  DELY  C  J ) 

CONTINUE 

RETURN 

END 


o  o  o 
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THIS  SUBPROGRAM  IS  USED  TO  PREDI CT  PRESSURE  AT  DISTANCE  X 
FROM  NOZZLE  ENTRANCE  (BY  EFFECTIVE  AREA)  IF  ISENTROPIC 
FLOW  HAS  BEEN  ASSUMED  ALONG  THE  NOZZLE 

SUBROUTINE  EXITP  <  R*  PO#  RD*  PEXIT  ) 

PE-0.001 

H-l. 

TOLER-0.00001 

DE-<  <R-  l>/2>**0.5*(  2/<R«-  !))**(  (R*  1)  /(  2*<R-  1)  )  >/RD 
RE-C l-(PE/PO)**<<R- 1)/R))**0.5*(PE/PG)**C 1/R) 

IF  ( ABS( < RE-DE) /DE) .LE. TOLER)  GO  TO  53 
IF  (RE.GT.DE)  GO  TO  1 
7  PEL-PE+H 

2  REL-<  I  -  (  PEL/PO) **<<R-l)/R)>-*0.5*<  PEL/PO) **( 1/R) 

IF  <ABS((REL-DE) /DE) .LE. TOLER)  GO  TO  53 

IF  ( (REL-DE) . GT.0)  GO  TO  3 
PE- PEL 
GO  TO  7 

3  H-H/2 
PEL- PE+H 
GO  TO  2 

1  PEL-PE+H 

10  REL-( l-<  PEL/PO)** <  <  R- 1 > /R) >**  0.5# ( PEL/PO)** ( 1/R) 

IF  < ABS( < REL-DE) /DE) .LE. TOLER)  GO  TO  53 
IF  <( REL-DE)  .GT.0)  GO  TO  5 
H-H/2 
PEL- PE*H 
GO  TO  10 
5  PE- PEL 

GO  TO  1 
53  CONTINUE 
PEXIT- PE 
RETURN 
END 


